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Abstract 



We used fuzzy cognitive mapping (FCM) to develop a generic shallow lake ecosystem model by augmenting the 
individual cognitive maps drawn by 8 scientists working in the area of shallow lake ecology. We calculated graph 
theoretical indices of the individual cognitive maps and the collective cognitive map produced by augmentation. 
There were a total of 32 variables with 113 connections in the collective cognitive map. The graph theoretical 
indices revealed internal cycles showing non-linear dynamics in the shallow lake ecosystem. The ecological processes 
were organized democratically without a top-down hierarchical structure. The most central variable in the collective 
map was submerged plants. The strongest connections were suspended solid concentration decreasing water clarity, 
phosphorus concentration increasing the phytoplankton biomass, higher water clarity increasing submerged plants, 
benthivorous fish biomass reducing submerged plants and increasing suspended solid concentration, and submerged 
plants decreasing suspended solids. 

The steady state condition of the generic model was a characteristic turbid shallow lake ecosystem. The generic 
shallow lake ecosystem model had the tendency to go into a turbid state since there were no dynamic environmental 
changes that could cause shifts between a turbid and a clearwater state, and the generic model indicated that only a 
dynamic disturbance regime could maintain the clearwater state. The model developed herein captured the empir- 
ical behavior of shallow lakes, and contained the basic model of the Alternative Stable States Theory. In addition, 
our model expanded the basic model by quantifying the relative effects of connections and by extending it with 22 
more variables and 99 more weighted causal connections. On our expanded model we ran 4 simulations: harvesting 
submerged plants, nutrient reduction, fish removal without nutrient reduction, and biomanipulation. Only biomanip- 
ulation, which included fish removal and nutrient reduction, had the potential to shift the turbid state into clearwater 
state. The structure and relationships in the generic model as well as the outcomes of the management simulations 
were supported by actual field studies in shallow lake ecosystems. Thus, fuzzy cognitive mapping methodology en- 
abled us to understand the complex structure of shallow lake ecosystems as a whole and obtain a valid generic model 
based on tacit knowledge of experts in the field. 
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1 Introduction 



Shallow lake ecosystems are one of the most numerous and one of the most vulnerable ecosystems in the world. Most 
shallow lakes are suffering from structural and functional changes and food web breakdown known as eutrophication 
as a result of excessive nutrient inputs from increased agricultural and urban activities, especially in highly popu- 
lated areas (Moss, 1998). Until the early 1990's, the management and restoration strategies for eutrophic shallow 
lakes focused on decreasing inflowing and in-lake nutrients, mainly phosphorus concentrations (Carpenter, Kitchell & 
Hodgson, 1985; Edmonson, 1985). However, eutrophic shallow lakes showed resilience and resistance to restoration 
(Sas, 1989; Scheffer et al, 1993a). The Alternative Stable States Theory showed that shallow lakes may be in two al- 
ternative stable states, a phytoplankton dominated turbid water state and a submerged plant dominated clearwater state 
(Scheffer et al. 1993a; Moss, 1990, 1998; Scheffer, 1998; Jeppesen, 1998). Each state is stable, in the sense that a shift 
from one state to the other requires a strong intervention from outside of the system. Thus, developing management 
plans for restoration and predicting the outcomes of those plans using mathematical modelling tools and computer 
simulations has been one of the key concerns of shallow lake ecologists and managers. So far, several models of 
shallow lake ecosystems have been developed. Many studies in the 1970s tried to model the eutrophication process in 
lakes using differential/difference equation models where parameters were empirically obtained (Ross, 1976, Lassiter 
and Kearns, 1977). These models assumed that the response of the system to the disturbance was in a linear fashion 
(Vollenweider, 1975), and were case-specific models. Later in the 1990's, more general models of shallow lakes were 
built, both aiming to explain the whole ecosystem structure (Scheffer et al., 1993a; Scheffer, 1998), and focusing on 
a particular compartment such as submerged plants (Scheffer, Bakema & Wortelboer, 1993b). The model developed 
by Scheffer et al. (1993a) became the basic model, which could be modified for specific cases. Recently, Zhang et al. 
(2003) employed a dynamic structural model for a shallow lake ecosystem. Although this latter model was successful 
in predicting the dynamics of the system, it was also case-specific. 

In this study, we use cognitive mapping to develop a shallow lake ecosystem model that is derived from the 
collective opinions and experiences of scientists who have been working on shallow lakes. The hypothesis is that this 
model will be a successful and detailed generic model of shallow lake ecosystems. We tested the results of this model 
with empirical results published in scientific literature. 

Cognitive maps are defined as qualitative models of how a given system operates. A cognitive map is based on 
defined variables, which can be physical quantities, or complex and abstract ideas, and the causal connections among 
them (Ozesmi and Ozesmi, 2004). Cognitive mapping is a useful tool to model complex systems such as shallow lake 
ecosystems. In addition, cognitive mapping allows for combining different expert opinions which result from different 
experiences in different physiogeographic contexts. 

2 Methods 

We collected the cognitive maps of 8 scientists and coded them as square matrices. To analyze the structure of 
the cognitive maps, we calculated the graph theoretical indices of the matrices. Then we augmented the matrices 
to obtain a collective cognitive map. The collective map was further used for simulations to test the outcomes of 
possible management options and the results were compared to the results of other studies in the literature to test the 
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performance of the model. 

2.1 Why use cognitive mapping? 

Cognitive mapping is a methodology where expert perception of complex systems can easily be obtained, combined 
and simulated. Cognitive mapping can model nonlinear relationships and feedbacks that are characteristics of many 
ecosystems. Although ecological applications are rare, the method has solid foundations. Cognitive maps are based 
on graph theory and digraphs (see Bang-Jensen and Gutin, 2002 for full treatment of the topic). Harary, Norman & 
Cartwright (1965) showed how directed graphs can be used to study the structural properties of the empirical world. In 
our study, we considered shallow lake ecosystems as empirical worlds. Since the 1970s, numerous studies have used 
cognitive mapping for representing complex systems and for decision-making (Axelrod, 1976; Bauer, 1975; Bougon 
et al., 1977, Brown, 1992; Carley and Palmquist, 1992; Cossette and Audet, 1992; Hart, 1977; Klein and Cooper, 
1982; Malone, 1975; Montazemi and Conrath, 1986; Nakamura, Iwai & Sawaragi, 1982; Roberts, 1973). In the 
1980s Kosko modified binary cognitive maps by applying fuzzy causal connections with real numbers in the range 
[—1,1] to the connections, using the term of fuzzy cognitive map (FCM), and provided means to simulate these FCMs 
(Kosko, 1986, 1987, 1988, 1992a, 1992b). The approach of FCM has been used in a wide variety of research, ranging 
from political developments (Taber, 1991) to data mining of the World Wide Web (Hong and Han, 2002; Lee et al., 
2002). Ozesmi (1999) used cognitive maps to represent local knowledge systems as told by informants in different 
stakeholder groups, the main assumption being that individuals have cognitive models that are internal representations 
of a partially observed world (Bauer, 1975). 

The use of FCM in ecology has been rather limited. Known examples include Radomski and Goeman (1996) who 
analyzed possible management options for sports fisheries; Hobbs et al. (2002), who used FCM to define management 
objectives for Lake Erie ecosystem; Ozesmi and Ozesmi (2003), who analyzed the perceptions of different stakeholder 
groups about a lake ecosystem in order to develop a management plan; and Skov and Svenning (2003), who combined 
FCM with a GIS to use expert knowledge to predict plant habitat suitability in a forest. 

2.2 Data Collection 

In this study, we used a multi-step fuzzy cognitive mapping approach (Ozesmi and Ozesmi, 2004; a detailed description 
of the standard method used can be found in this methodological paper) to model the system dynamics of shallow 
lakes based on models drawn/mapped by 8 expert scientists in the discipline. The scientists were all male, with an 
average age of 37.9 ± 7.3 years. The question they were asked was 'What are the variables, factors determining 
the state of shallow lake ecosystems and how do they interact with each other?' They defined variables and drew 
causal relationships among these variables. In addition, they indicated the relative strength/weight and direction of the 
relationships using a scale varying between -1 and +1 with 0.25 increments. The sample size (i.e., the number of maps) 
was considered sufficient once the number of new variables leveled off as the maps were added together (Ozesmi and 
Ozesmi 2004, p. 48). The number of new variables added per new map leveled off at about 1 by the 8 th interview after 
1000 Monte-Carlo simulations (figure 1) using Estimates 6. Obi (Colwell, 1997). Accumulation of variables indicated 
that 8 maps were sufficient to characterize the system under study. 
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2.3 Data Analysis 



Once the maps were drawn, they were coded as square (adjacency) matrices according to graph theory in the form 
A(D) = [a,ij] (Harary et al., 1965). The variables v t are listed on the vertical axis and Vj on the horizontal axis. When 
a connection exists between two variables the value is coded in the square matrix (between -1 and 1). 

Then the structure of the maps was analyzed using indices derived from graph theory. We calculated the indegrees 
and outdegrees of each variable for each individual map. Indegree is defined as the column sum of absolute values of 
a variable and shows the cumulative strength of variables entering the variable. Outdegree is defined as the row sum 
of absolute values of a variable and shows the cumulative strengths of connections exiting the variable. In addition, 
we determined the variables with positive outdegree and zero indegree, with positive indegree and zero outdegree and 
with non-zero indegree and outdegree, which represent transmitter, receiver and ordinary variables, respectively. For 
each individual map to determine its network structure, we also calculated the centrality, which is the sum of indegree 
and outdegree and indicate the contribution of that variable to the cognitive map; complexity index, which is the ratio 
of number of receiver variables to that of the transmitter variables; connectivity as the number of connections; density, 
which is calculated as a function of connectivity, and shows how connected or sparse the maps are; distance, which 
is the average causal link between variables and hierarchy index, which is a function of the outdegrees and number 
of variables in a given map and shows how democratic the map is, (see Ozesmi and Ozesmi, 2004, p. 50-51. for 
formulas). 

Individual maps were later augmented (matrices were added together in an overlapping manner (Kosko, 1987, 
1992a,b)) to obtain a collective cognitive map (social map) on the dynamics of shallow lake ecosystems. First an 
augmented matrix is created that includes all the variables from all the individual cognitive maps. Second each 
individual cognitive map is coded into an augmented matrix and the individual cognitive maps are added together 
using matrix addition to create a social cognitive map. As the maps are added together, conflicting connections 
with opposite signs will decrease the causal relationship, while agreement reinforces causal relationships, forming a 
collective map, which is considered a consensus opinion of the scientists. While augmenting, equal weight was given 
to the map of each scientist. 

The collective cognitive map was used in simulating system behavior and to run management simulations. Sim- 
ulations were made by multiplying a vector of initial states of variables (I n ) with the adjacency (square) matrix A of 
the collective map. Usually a threshold function or a transformation by a bounded monotonic increasing function is 
applied to the result of the matrix multiplication, I n x A, at each simulation time step (Kosko 1987; Kosko 1992b, 
Ozesmi and Ozesmi 2004, p. 55). We used a logistic function (1 + e^) -1 to transform the results into the interval 
[0, 1]. The resulting transformed vector was then repeatedly multiplied by the adjacency matrix and transformed until 
the system converged to a fixed point. We produced 10 random initial state vectors and ran the model to determine the 
stability of the collective map without any management. Later management options were simulated by clamping given 
variables through the multiplication process. Relevant variables were selected based on common practice shallow lake 
management interventions and these were simulated. The simulation results were compared with results of empirical 
studies. All computations were carried out using subroutines implemented in Octave programming language version 
2.0. 16 on a Linux OS. 
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3 Results 



The mean number of variables in the individual cognitive maps of the shallow lake ecosystems drawn by 8 scientists 
was 14.6 ± 3.3 SD, while there were 20.9 ± 5.7 connections on average between the variables that they defined. 
There were a total of 32 variables with 113 connections in the collective cognitive map obtained by augmenting the 8 
individual maps (table 1). 

Mean values of graph theoretical indices of the individual maps and the collective map are shown in table 1. 
Scientists have 3.8 ± 2.5 SD transmitter variables in their individual maps, whereas in their collective cognitive map, 
there are a total of 6 transmitter variables. These were external nutrient load, micropollutants, dissolved organic carbon 
concentration, temperature, fetch and shoreline slope. There were about the same number of receiver variables (utility 
variables), 3.4 ±2.1 SD in the individual maps, whereas only one in the collective map (sediment phosphorus release). 
There were 10.6 ±3.1 ordinary variables (means) in the individual maps, and 25 in the collective map. 

Looking at the cognitive maps, we observe that there are on average 1.56 ± 0.44 SD causal chains in the individual 
maps, whereas there are 2.79 causal connections in the collective map. The complexity, density and hierarchy indices 
of the cognitive maps are low (table 1). 

The most frequently mentioned variables that are recurrent in the 8 cognitive maps are given in table 2. Suspended 
solids were mentioned in each map, and water clarity, zooplanktivorous fish, submerged plants, benthivorous fish and 
phosphorus concentration were mentioned in 7 of the 8 maps. 

The most central variable in the collective map was submerged plants (figure 2). Moreover, submerged plants have 
as much effect on other variables (outdegree: 2.8) as it is affected by others (indegree: 3.0). 

In figure 2, variables are ordered according to their centralities. Variables shown in the figure have centrality 
higher than 1 .0, and were mentioned by more than half of the participating scientists. Suspended solid concentrations, 
phytoplankton biomass, water clarity and zooplankton biomass were more affected by other variables than affecting 
others (indegree > outdegree), while the situation was the opposite for phosphorus concentration, benthivorous and 
planktivorous fish biomass, waves and wind action, benthic biomass and water level. 

We drew the cognitive interpretation diagram (CID, figure 3) according to Ozesmi and Ozesmi (2004, p. 53). 
This diagram transforms the complex "spaghetti and meat balls" graph of the collective map into an understandable 
format where connections with weights more than 0.25 are drawn, scaled to their actual weights. Onto our CID, we 
superimposed the basic model of Scheffer et al. (1993a). Variables mentioned also by Scheffer et al., are drawn in 
two concentric circles, compared to those with only one circle, which are our additional variables. Binary connections 
that are also present in Scheffer et al. (1993a) have a circle at the arrowheads and have all the same signs as ours. 
Connections missing in our model (from submerged plants to allelophatic substances, and from there to phytoplankton, 
from submerged plants to zooplankton-fish interaction) have a thin line with a circle at the pointing end. 

Accordint to the CID (figure 3), the strongest connections are suspended solid concentration decreasing water 
clarity, phosphorus concentration increasing the phytoplankton biomass and higher water clarity increasing submerged 
plants (see table 3 for actual weights). These connections are followed in strength by benthivorous fish biomass 
reducing submerged plants and increasing suspended solid concentration, submerged plants decreasing suspended 
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solids and zooplankton density reducing the phytoplankton biomass (table 3). 

We determined the steady state condition of the system before any management option was considered, in order to 
be able to see the tendency of the ecosystem based on the collective cognitive map. We ran the generic shallow lakes 
ecosystem model with 10 different random initial states for all variables between and 1 from a uniform distribution. 
We also ran an initial state vector where all variables were set to 1, as well as with a real data set compiled from Lake 
Eymir database (Burnak and Beklioglu, 2000; Tan, 2002). In all of these no management simulations, the system 
reached steady state after 10 - 12 iterations, where the states were the same to 10 -6 digits. A one-way ANOVA test 
did not show any significant differences (Fn.341 : 0.000, p : 1.000). The mentioned steady state condition of the 
system is given in figure 4. The steady state condition of the collective generic model is a characteristic turbid shallow 
lake ecosystem where primary production was governed by phytoplankton and benthic primary production. As figure 
4 clearly shows, phytoplankton biomass, phosphorus, and suspended solids concentrations are the highest variables, 
and therefore not surprisingly, water transparency and submerged plants are almost 0. 

Following the determination of system's steady state, management scenarios are developed based on common 
interventions in the field of shallow lake restoration. These scenarios were: A) Harvesting submerged plants; B) 
Nutrient reduction; C) Fish removal without nutrient reduction; and D) Biomanipulation involving both fish removal 
and nutrient reduction. The results of the managemnt simulations are shown in figure 5. We looked how much the most 
central variables (water clarity, phosphorus concentration, phytoplankton biomass, zooplankton density, planktivorous 
fish biomass, benthivorous fish biomass, benthic biomass, submerged plants, suspended solid concentration, water, 
level, and waves/wind action) changed from their steady state values. A negative value indicates a reduction in the 
variable state compared to its steady state, while a positive value reflects an increase in the variable state. 

Harvesting Submerged Plants. Submerged plants are usually harvested in recreational lakes where it is considered 
a nuisance to boating and angling as well as aesthetically displeasing. When submerged plants are harvested (figure 5- 
A), the main change observed is an increase in the suspended solid concentration. While the phosphorus concentration 
and phytoplankton biomass slightly increased, zooplankton density decreased. An interesting result is the increase 
in the wave and wind action, revealing the importance of submerged plants for stabilizing the wave action, and thus, 
decreasing the suspended solids. When we look at the scale of change in variables under this scenario, compared to 
other management options (figure 5-B, C, D), we see that it is relatively low. The low relative change is reasonable, 
considering that the system at its steady state is already in turbid water conditions with very low water transparency 
and few submerged plants (figure 4). 

Nutrient Reduction. In this simulation (figure 5-B), nutrient reduction leads to a strong decrease in the phytoplank- 
ton biomass, a reduction in the concentration of suspended solids, and a decrease in the planktivorous and benthivorous 
fish biomass. 

Fish Removal without Nutrient Reduction. When the zooplanktivorous fish were removed, naturally the zooplank- 
ton density increased and the phytoplankton biomass decreased. By also removing benthivorous fish, the suspended 
solids decreased considerably in our model (figure 5-C). 

Bimanipulation. In the last simulation (figure 5-D) benthivorous and planktivorous fish were removed and nutrients 
were also reduced. In comparison to all the other scenarios we observed a much higher change in the variables and a 
shift towards a clearwater state. For the first time submerged plants increased considerably as well as water clarity and 
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zooplankton density, while phytoplankton biomass decreased. Suspended solids also decreased mainly as an outcome 
of the decreased benthivorous fish stock, as suggested in scenario 3 (fish removal). 

4 Discussion 

4.1 Graph Theoretical Indices 

The low mean number of transmitter variables in the cognitive maps of scientists indicates that they define relatively 
few outside controls (forcing functions). The mean number of receiver variables, being also low, and the high number 
of ordinary variables in the cognitve maps further show that the shallow lake ecosystem is predominantly governed by 
circular internal ecological processes. These internal processes define how the ecosystem will behave. Furthermore, 
the low values of the complexity, density and hierarchy indices, considered among with the mean number of causal 
connections, indicates that according to the scientists, the processes in the shallow lake ecosystem are organized 
"democratically" without showing a top-down hierarchical structure, and there are many internal cycles showing non- 
linear dynamics. 

4.2 Variable characteristics 

Although the most mentioned variables show what is in the minds of the participating scientists, we can understand 
the most important variables that control the dynamics of the shallow lake ecosystem by looking at the centrality of 
variables in the collective cognitive map. The indegree and outdegree, which add up to centrality, will indicate how 
much the variable is affected by other variables or how much it affects the others. Thus, the differences in indegrees 
and outdegrees of the variables are important for inferring the relative role of that variable in the ecosystem. 

The most central variable in the collective map is found to be the submerged plants, assigning a key central role to 
submerged plants. This result corroborates the field studies and alternative stable states theory, where submerged plants 
are the key components in determining the ecosystem dynamics in a shallow lake (Scheffer et al., 1993a; Jeppesen, 
1998). Moreover, close values of indegree and outdegree of submerged plants indicated that they have as much effect 
on other variables as they are affected by others. This result is similar to the figure of the basic shallow lake ecosystem 
model by Scheffer et al. (1993a), where vegetation has a central position and both affects and is affected by other 
variables. 

Variable centralities in the cognitive map (figure 2) revealed two dominant groups of variables. Accordingly, 
suspended solid concentrations, phytoplankton biomass, water clarity and zooplankton biomass were more affected 
by other variabes than affecting others, whereas the oppoite is true for phosphorus concentration, benthivoros and 
planktivorous fish biomass, waves and wind action, benthic biomass and water level. The latter group will act more like 
forcing functions, influencing ecosystem processes (high outdegrees) while being less influenced by other ecosystem 
components (low indegrees). While these results are not surprising according to our current knowledge of shallow lake 
ecosystems, it is surprising that the generic model captures the relationship of the benthi-planktivorous fish biomass 
with other components of the ecosystem. As the model shows, fish will be less affected by other variables having a 
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much longer life span compared to other organisms, such as zooplankton. During their longer life span, fish will have 
a lasting effect on the shallow lake ecosystem dynamics. 

4.3 Cognitive Interpretation Diagram of the Shallow Lake Ecosystem 

Cognitive interpretation diagram (CID, figure 3) revealed the fact that basic shallow lake model introduced by Scheffer 
et al. (1993a) is consistent with and included in our model. In our model, allelophatic substances do not have a strong 
effect in the system. The advantage of our model is that all connections regardless whether they are included in 
Scheffer et al. (1993a) or not, have weights associated with them showing their relative importance in the ecosystem 
dynamics. 

The strength and pattern of the connections in our model, according to CID, are also consistent with actual field 
studies in shallow lake ecosystems and theoretical conclusions derived from these (Jeppesen, 1998; Moss, 1998; 
Scheffer, 1998). Our model contains the basic model of Scheffer et al. (1993a). This is not surprising because all of the 
scientists that contributed their cognitive maps certainly knew of Scheffer's model and were influenced by it. However 
the collective cognitive map improves Scheffer's basic model by quantifying the relative effects of connections and 
adding 22 more variables and 99 causal connections. 

4.4 Simulations Using Collective Fuzzy Cognitive Map 

The multi-step fuzzy cognitive mapping approach (Ozesmi and Ozesmi, 2004) enables us to run the collective cog- 
nitive map as a generic shallow lake ecosystem simulation model. We can model the system as it is without any 
management, run simulations of management options by fixing state values of relevant variables and find the impact 
of those management options by looking at deviations from the steady state. One should be careful in interpreting 
results obtained from simulations of fuzzy cognitive maps, because these are not system dynamics models in the 
classical sense, yielding directly measurable results. Rather, the results are qualitative indicators showing whether a 
variable increases or decreases relative to its steady state or to other variables in the system. 

4.4.1 Steady State Condition 

Although we did not indicate anything regarding the state of the system to the 8 scientists who drew their maps, their 
collective generic model, where there is no human management, ends up at a steady state condition that is in turbid 
water state. Moreover, we did a simulation to test this steady state, excluding all the possible variables that may be 
considered as outside disturbances to the system (external nutrient load, sediment phosphorus release, micropollutants, 
bioturbation, and waves and wind action), and the system still stabilized on a turbid water state. This might be because 
these scientists are experts on management and restoration of shallow lake ecosystems which are usually already in 
turbid water state. Alternatively, this result raises the question whether the natural tendency of most shallow lake 
ecosystems is to stay in a turbid state if conditions are stable over a long period of time. Maybe, what maintains or 
causes shifts from turbid to clearwater state is a dynamic environment in which variables change based on conditions 
outside the system. In other words, is it a dynamic disturbance regime that maintains the clearwater state in shallow 
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lake ecosystems? We know that Fox (1979), Connell (1979), and Reynolds, Padisak & Sommer (1993) postulate in the 
intermediate disturbance hypothesis that high species diversity in an assemblage is maintained by fluctuations in the 
environmental conditions occurring within the generation time of those organisms forming the assemble. Similarly, 
our model suggests that diversity at the ecosystem level is supported by outside disturbances. This, in turn, would 
contribute to diversity at the species level. 

Our shallow lake ecosystem model represents the ecosystem processes as they are. Thus applying different man- 
agement options can be considered as disturbing the system in a directed manner, taking advantage of the system 
characteristics to maintain a clearwater state. Maintaining a clearwater state is preferred for reasons of productivity, 
biodiversity and aesthetic concerns. But if all shallow lakes in the world were maintained in a clearwater state this 
would be a loss in total biodiversity. However, today, human economic activities cause most shallow lakes to turn into 
turbid ecosystems. Therefore, we also explored different management options that can provide for a clearwater state. 

4.4.2 Simulation of Management Scenarios for Shallow Lakes 

The collective cognitive map of the shallow lake ecosystem can be considered an a priori model, because no empirical 
studies were used to construct or parameterize the model. Of course, one could argue that the scientists who contributed 
to the generic model have done so based on their own empirical work and knowledge derived from the existing 
literature. However, the models that they shared with us are based on tacit knowledge that does not involve explicitly 
measurable relationships but qualitative judgments on cause and effect relationships between variables. By comparing 
the performance of this a priori qualitative model directly with empirical results, we can continue to test its generic 
validity. 

Harvesting Submerged Plants. This scenario agrees with the general conclusions of field studies where the loss 
of submerged plants was followed by an increase in suspended solid concentrations mainly because of an increase in 
wind and wave action (Bostrom, Jansson & Forsberg, 1982). Moreover, Breukelaar et al.(1994) reported that the loss 
of submerged plants leads to an increase in the benthic fish biomass. They concluded that it is easier for benthivorous 
fish to feed when there are no submerged plants. This is exactly what we have observed in our simulation (figure 5- A). 
Moreover, it is known that the submerged plants enhance the zooplankton density by providing refuge to zooplankton 
against zooplanktivorous fish predation (Timms and Moss, 1984). As submerged plants have been removed in our 
simulation, we observe a drop in zooplankton density. The expected increase in phytoplankton biomass due to a 
reduction in grazing zooplankton density is also observed in this simulation. 

Nutrient Reduction. The strong decrease in the phytoplankton biomass and reduction in the concentration of 
suspended solids is probably as a result of the observed decrease in the planktivorous and benthivorous fish biomass. 
In agreement with alternative stable states theory (Scheffer et al., 1993a), despite the strong decrease in phytoplankton 
biomass, we do not observe any improvement in the water clarity and submerged plants which can be attributed to the 
resilience of the system. Field studies have corroborated this resilience repeatedly. (Scheffer, 1998; Breukelaar et al., 
1994; Jeppesen et al., 1997, 1999; Jeppesen, 1998; Beklioglu, Carvalho& Moss, 1999). 

Fish Removal without Nutrient Reduction. Our model is in agreement with Scheffer (1998). When th zooplanktiv- 
orous firs are removed, naturally the zooplankton density increases, and as a consequence, the phytoplankton biomass 
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decreases due to intensified grazing (Scheffer, 1998). Also in our model, suspended solids were decreased by remov- 
ing benthivorous fish species. Benthivorous fish species are known to enhance suspended solids concentration due 
to their feeding action on the sediment (Meijer et al., 1990; Breukelaar et al., 1994; Sidorkewicj, Lopez Cazorla & 
Fernandez, 1996). When benthivorous fish feed on the sediment, they stir up the phosphorus bound in the sediment 
into the water column. (Horpilla and Kairesalo, 1990; Breukelaar et al., 1994). We observed this effect by a slight 
decrease in the concentration of phosphorus when benthivorous fish were removed. 

We did not observe,however, an increase in water clarity despite the reduction in suspended solids and phytoplank- 
ton biomass when only fish were removed. The same has been observed in empirical studies (Jeppesen et al, 1997, 
1999; Hansson et al., 1998; Beklioglu et al., 1999). The alternative stable states theory (Scheffer et al., 1993a) explains 
this phenomenon as the resilience of the system to shift into a clearwater state. The resilience of the turbid water state 
is a result of the delayed development of submerged plants. Before submerged plants can get truly established in the 
lake and their stabilizing buffer mechanisms can take effect, high nutrient concentrations in the water favor phyto- 
plankton redevelopment (James & Barko, 1990; Scheffer et al., 1993a; Weisner, Strand & Sandsten, 1997; van Donk 
& van de Bund, 2002). 

Biomanipulation. The shift to clearwater state in this scenario is consistent with the biomanipulation literature 
(Shapiro, Lamarra & Lynch, 1975; Benndorf et al., 1988; Hosper, 1989; Jeppesen et al., 1990; Carpenter and Kitchell, 
1993; Philips and Moss, 1994; Hansson et al., 1998; Meijer, 2000). Moreover, the action of waves and wind decreased, 
most probably as an outcome of submerged plant development (Bostrom et al., 1982). 

In summary, we used cognitive mapping to develop a generic shallow lake ecosystem model that is derived from 
the collective tacit knowledge and experiences of 8 scientists who have been working on shallow lakes. There were a 
total of 32 variables with 113 connections in the collective cognitive map. The graph theoretical indices show that the 
shallow lake ecosystem is predominantly governed by internal cycles showing non-linear dynamics. The ecological 
processes are organized democratically without showing a top-down hierarchical structure. 

Submerged plants were the most central variable in the collective map. The strongest connections in decreasing 
order were: suspended solid concentration decreasing water clarity; phosphorus concentration increasing the phyto- 
plankton biomass; water clarity improving submerged plants; benthivorous fish biomass reducing submerged plants 
and increasing suspended solid concentration; submerged plants decreasing suspended solids; and zooplankton den- 
sity reducing the phytoplankton biomass. The steady state condition of the generic model was a characteristic turbid 
shallow lake ecosystem where primary production was governed by phytoplankton and benthic primary producers. 
The turbid state might be because participating scientists are experts on management and restoration of shallow lake 
ecosystems that generally show a high trophic state. Alternatively, most shallow lake ecosystems have the natural 
tendency to stay in a turbid state if conditions are stable over a long period of time. Shifts between a turbid and a 
clearwater state may be the result of a dynamic environment in which variables change in relation to the conditions 
outside the system. Our model indicates that only a dynamic disturbance regime could maintain the clearwater state in 
shallow lake ecosystems similarly to how species diversity is maintained in the intermediate disturbance hypothesis. 
Diversity at the ecosystem level created by outside disturbances can contribute to diversity at the species level in large 
physiogeographic scales. 

The generic FCM model captured all the behavior predicted by the Alternative Stable States Theory developed for 
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shallow lakes ecosystems. Our model contains the basic model of Scheffer et al. (1993a) but expands it by quantifying 
the relative effects of connections and adds 22 more variables and 99 more weighted causal connections. On our 
expanded model we ran 4 simulations: harvesting submerged plants, nutrient reduction, fish removal without nutrient 
reduction, and biomanipulation. Harvesting submerged macrophytes did not show a strong impact on the system 
because the steady state of the system was already turbid. This simulation showed that submerged plants have an 
important effect even in a turbid state where they are minimal. Similarly, neither nutrient reduction nor fish removal 
was enough by itself for the lake to shift toward a clearwater state. However, fish removal accompanied with nutrient 
reduction (biomanipulation scenario) showed that this strategy is applicable for restoring turbid water systems. Not 
only the structure and relations in the generic model were supported by actual field studies in shallow lake ecosystems 
and theoretical conclusions derived from these (Jeppesen, 1998; Moss, 1998; Scheffer, 1998) but also all the outcomes 
of the 4 management simulations were supported by other empirical results in literature. 

Cognitive mapping methodology also enabled us to explore a complex system in detail and gain insights about how 
the ecosystem functions without being trapped into the specificity of sites but be able to generalize from them. The 
model produced outcomes consistent with empirical results. The behavior of the model suggested that shallow lake 
ecosystems have a tendency to become turbid unless there are dynamic forcing functions changing central variables. 
This theoretical insight is consistent with the theory that there is a tendency for an ecosystem to deteriorate unless the 
dynamics of outside forces rejuvenate critical central variables (Connell, 1979; Fox, 1979; Reynolds et al., 1993). 
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Table 1. Average (± SD) graph theoretical indices of the individual cognitive maps and the indices of the collective 
cognitive map of the scientists. 
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Table 2. Most mentioned variables listed in decreasing order of number of times occurring in the maps drawn by more 
than two scientists. 
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Table 3. The strongest connections in the collective cognitive map. The relations are given in decreasing weight. The 
direction of the connection is from variable 1 to variable 2. 
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Figure Captions 



Figure 1. Accumulation curve of number of new variables versus the numer of added maps. The dashed lines represent 
±SD. 

Figure 2. Most central variables in the collective map. Bars indicate centrality, which is the sum of the indegree and 
outdegree of a variable showing the importance of the role of the variable. 

Figure 3. The cognitive interpretation diagramm (CID) of the collective map of 8 scientists. Only the strongest 
connections (> 0.25) are drawn. The thickness of the lines indicates the relative strength of the connections (2 pixel 
corresponds to a connection with strength 0.125). Solid lines represent positive connections, dashed lines negative 
connections. The generic model of Scheffer et al. (1993a) is superimposed on our CID and is represented by two 
concentric circles. Binary connections that are also present in Scheffer et al. have a circle at the arrowheads and have 
all the same signs as ours. 

Figure 4. imulation of steady state conditions of the variables obtained from the collective cognitive map with no 
management option. 

Figure 5. The change from steady state under different management scenarios. Scenario A: Submerged plant harvest; 
B: Nutrient reduction; C: Fish removal without nutrient reduction; D: Biomanipulation, fish removal accompanied 
with nutrient reduction. 
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